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"^ On large scales, the higher order moments of the mass distribution, Sj = 0/^2 i ^-S-i the skewness 

■^ 5^3 and kurtosis S4, can be predicted using non-linear perturbation theory. Comparison of these pre- 

O , dictions with moments of the observed galaxy distribution probes the bias between galaxies and mass. 

Applying this method to models with initially Gaussian fluctuations and power spectra P{k) similar 
to that of galaxies in the APM survey, we find that the predicted higher order moments Sj{R) are in 
good agreement with those directly inferred from the APM survey in the absence of bias. We use this 
result to place limits on the linear and non-linear bias parameters. Models in which the extra power 
observed on large scales (with respect to standard CDM) is produced by scale-dependent bias match 
the APM higher order amplitudes only if non-linear bias (rather than non-linear gravity) generates the 
observed higher order moments. When normalized to COBE DMR, these models are significantly ruled 
out by the ^3 observations. The cold plus hot dark matter model normalized to COBE can reproduce 
the APM higher order correlations if one introduces non-linear bias terms, while the low-density CDM 
model with a cosmological constant does not require any bias to fit the large-scale amplitudes. 

Subject Headings: Large-scale structure of the universe — galaxies: clustering 



1 Introduction 

In standard models for galaxy formation, such as cold dark matter (CDM) and its offshoots, it is 
usually assumed that the observable galaxy distribution is related through a simple bias mechanism to 
the underlying matter distribution predicted by theory (e.g., Bardeen,e< al. 1986). In essence, following 
Kaiser (1984a, b) and Bardeen (1984), one assumes that galaxies form from peaks above some global 
threshold in the smoothed linear density field. In the limit of high threshold and small variance, the 
peaks model is well approximated by the commonly employed linear bias scheme, in which the galaxy 
and mass density fields, Sg(x) = (ng(x) — ng)/ng and S(x) = {p{x) — p)/p, are related through a 
constant bias factor, 

Sg(x) = hgS{x) . (1) 

More generally, if the bias is local^ the galaxy field may be some non-linear function of the density field, 
Sg{x) = f{6{x)). Recently, more complex, non-local schemes for biasing have also been studied (Babul 
and White 1991, Bower, et al. 1993); in these models, the effective bias factor becomes scale-dependent. 
If the bias factor increases with scale, the galaxy spectrum will have more relative power at large scales. 
Indeed, it has been shown that such modifications of the standard bias scheme can generate the excess 
large-scale power found, e.g., in the APM survey (Maddox et al. 1990) within the context of the standard 
(Qh = 0.5) CDM model. This modification of the standard CDM scenario is fundamentally different 
from other CDM variants, such as non-zero A, tilt, or cold plus hot dark matter, where there is genuine 
extra power in the density field. 

To date, two-point statistics (correlation functions and power spectra) have been the dominant 
benchmark for testing theories of biased structure formation. It is difficult, however, to test biasing 
models with two-point statistics alone, because one must specify a number of input variables to predict 
the amplitude of the power spectrum, e.g., the density of the universe Q, the Hubble constant h (where 
Hq = IQO/i km s~^ Mpc~^), the type of dark matter, the cosmological constant A, the initial ampli- 
tude and the clustering interactions (gravity). On the other hand, the non-linear perturbation theory 
predictions for the higher order moments of the mass distribution, Sj = 0/^2 i ^^^ approximately 
independent of time, scale, density, or geometry of the cosmological model (Juskiewicz, Bouchet, and 
Colombi 1993, Bernardeau 1994). The higher moments only depend on the hypothesis that the initial 
fluctuations are small and quasi-Gaussian and that they grow via gravitational clustering. Thus, the 
moments Sj are an excellent tool for discriminating among biasing scenarios. 

Frieman & Gaztahaga (1994) have considered how the higher order irreducible moments of the 
galaxy distribution can be used as a test of models for large-scale structure. In particular, they studied 
the three-point function in the standard CDM model and its variants with extra large-scale power 
(e.g., Clh = 0.2 CDM), as well as in a generalized version of the non-local, scale-dependent bias scheme 
embodied in the cooperative galaxy formation (hereafter CGF) model of Bower et al. . Using the results 
of second-order perturbation theory (Fry 1984), they compared the predictions of these models for the 
three-point function ^3 with data from the Center for Astrophysics (Huchra, et al. 1983), Southern Sky 



(Da Costa, e< al. 1991), and Perseus-Pisces (Haynes and Giovanelli 1988) redshift surveys in the mildly 
non-linear regime (^2 < !)■ 

In this Letter, we compare these model predictions to the higher order galaxy correlations recently 
estimated for the APM survey by Gaztahaga (1994). The values for the Sj moments in the APM survey 
are statistically more reliable than those inferred from the redshift surveys above, because the APM 
catalog covers over 200 times more volume and contains 300 times more galaxies than the combined 
CfA/SSRS. The skewness ^3 found in the APM survey is typically larger than in the optical redshift 
surveys (Gaztafiaga 1992). However, the disagreement is only significant on small scales: at the largest 
scales probed by the redshift surveys, R ~ 20 h~^ Mpc, the value of ^3 -^ 2 found in the CfA/SSRS 
catalogs is in good agreement with the APM results. Furthermore, comparison of the CfA and SSRS 
catalogs with comparable sub-samples of the APM indicates that the difference in skewness on small 
scales is attributable to the fact that the local galaxy distribution traced by the CfA and SSRS catalogs 
does not correspond to a "fair sample" rather than to uncertainties in the APM selection function 
(Gaztafiaga 1994). 

2 The APM galaxy moments 

In the weakly non-linear regime (^2 < l)i non- linear perturbation theory has been shown to hold remark- 
ably well when compared with N-body simulations (Baugh, Gaztafiaga, & Efstathiou in preparation, 
Efstathiou et al. 1988, Bouchet and Hernquist 1992, Weinberg & Cole 1992, Juskiewicz, Bouchet, & 
Colombi 1993, Bernardeau 1994). In perturbation theory (Fry 1984, Goroff, et al. 1986, Bernardeau 
1992, 1994), the matter density field evolved gravitationally from Gaussian initial conditions leads to 
a hierarchical clustering pattern of the form ^j = Sj ^2 ■ Here ^j is the volume-averaged J— point 
correlation function, 

lj{V) = ^ j ■■■ j d\i...d\jij{n,...,Tj)W{ri)...W{rj) , (2) 

where W{r) is a window function of characteristic volume V (below, we use a top-hat window of radius 
R). For models without strong features in the linear power spectrum (e.g., CDM and its popular 
variants), the perturbative J— point moments Sj depend only weakly on the window smoothing scale 
(for power law spectra, the Sj{R) are independent of R). Here we focus on the hierarchical amplitudes 
5^3 and ^4, the normalized measures of the skewness and kurtosis of the smoothed density probability 
distribution function. 

Figure 1 shows the values of S'i{R) and S4{R) for spherical cells of radius i?, estimated from the 
angular amplitudes s^{d) and 84(0) in the APM galaxy survey (cf. Fig. 4 of Gaztafiaga 1994). The 
error-bars correspond to the one-sigma dispersion between four disjoint zones of the APM map and 
are therefore an estimate of the 'cosmic variance'. The bins in the data correspond to the bins in the 
angular data: the original angular cell counts are done in conical volumes. Each angular bin in the 
sj{9) data is such that the conical volume in one cell is at least 100% larger than the volume of the 
next smallest bin, so the correlations in each bin are roughly independent from the others. The angular 
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Figure 1: The volume-averaged skewness S3{R) and kurtosis S4{R) for spherical cells of radius R = OD, 
estimated from sj(6) in the APM survey, where the effective depth V ~ 400 /i~^ Mpc. Error-bars 
correspond to one-sigma dispersion of sj{9) between 4 zones of the APM map. Curves correspond 
to the Sj{R) predictions in perturbation theory for standard CDM (T = 0.5, long-dash curves), low- 
density CDM (r = 0.2, short-dash curves), the inferred APM spectrum (solid curves) and standard 
CDM modified by the scale-dependent bias (CGF model), with k = 2.29, Rs = 20 /i"^ Mpc (dotted 
curve). These models are all shown with &i = 6 = 1, C2 = C3 = 0. Also shown (dot-dash curves) are the 
moments for a biased P = 0.2 model with b = 1.5, C2 = 0.45, C3 = 0.2 choosen to give similar predictions 
to the unbiased case. 



amplitudes have been corrected by a factor 0.95 to account for residual contamination from merged 

images in the APM survey (Maddox, et al. 1990, Gaztafiaga 1994); since this factor is approximate, 

the error-bars on sj{0) were constrained to be at least as large as the correction factor. We have also 

performed N-body simulations to verify that the errors introduced by uncertainties in the APM selection 

function and in the two- to three-dimensional inversion factors are small compared with those above 

(Gaztafiaga 1994, and in preparation). Thus, for the purposes of comparison with theoretical models, 

the data points shown in Fig. 1 can be taken to be independent and have been assigned conservative 

errors. (For additional details, see Gaztafiaga 1994.) 

For the models, we assume the initial fluctuations are small and Gaussian, and that they evolve under 

the influence of gravity in an expanding universe. We estimate the moments Sj(R) = ^j(R)/^2{^) 

for the matter distribution using (J — l)-order perturbation theory (Fry 1984, Bernardeau 1994). The 

results for Sj are completely determined by the linear power spectrum P{k), for which we adopt the 

parametric form 

2, A . 1.7k _ 9fc3/2 fc2\-2 



P(k) = Aaikil + ^^ + ^;^ + ^ 



(3) 



where the wavenumber k is in units of h Mpc^ , and ag = £,2i-^ ~ 8 h^ Mpc) is the variance of the 
linear mass fluctuation within (top-hat) spheres of radius 8 /i^ Mpc. For CDM models, equation (3) 
gives an accurate fit to the linear power spectrum with P = Qh (Davis, etal. 1985). We consider two 
models: 'standard CDM' with P = 0.5, and 'low-density CDM' with P = 0.2. For each model, we use the 



power spectrum P{k) to calculate the variance £,2i^) ^^"^ its derivatives 7„ = d" log ^2(-^)/'^log" R s-nd 
use the results of Bernardeau (1994) to calculate the predicted amplitudes S3{R) and S4{R). These are 
shown as the long and short-dash curves in Fig. 1 for standard and low-density CDM. For comparison, 
the solid curves in Fig. 1 show the predicted higher moments using the power spectrum inferred directly 
from the APM catalog by Baugh & Efstathiou (1993) as the linear input for higher order perturbation 
theory. This spectrum is reasonably well approximated by the T = 0.2 CDM model. 

To show the effects of scale-dependent bias, we also consider the non-local bias model of Bower, etal. 
(1993). In this model, the linear galaxy power spectrum is just the standard T = 0.5 CDM spectrum 
multiplied by the 'cooperative' bias factor b'^(k) = (1 + Ke~^ '' '^)^. In the parameter range studied 
by Bower, etal., the choice k = 2.29, Rg = 20 /i~^ Mpc for the strength and range of cooperative effects 
gives the best fit to the APM angular correlation function, and for illustration we focus on this case. 
The corresponding skewness and kurtosis are shown as the dotted curves in Fig. 1. As noted by Frieman 
& Gaztanaga (1994), in this model there is a sharp downturn in the Sj on large scales, i? > 10 /i~^ Mpc, 
which is at variance with the relative flatness of the observed moments. 

The model curves in Fig. 1 were all plotted assuming the linear bias model of equation (1) with 
6=1. When considering higher-order perturbations, however, one should self-consistently allow for 
higher-order (non-linear) bias, and replace (1) with an expansion of the form 

k\ 



^9 = fi^) = Y.T/ ■ (4) 



Fry and Gaztanaga (1993) have shown that such a local transformation S ^ 6g preserves the hierarchical 
nature of the matter distribution in the limit of small fluctuations, i.e., on large scales, although the 
values of the resulting galaxy hierarchical amplitudes Sj will in general differ from those of the density 
field (see also Juszkiewicz, etal. 1993). For J = 3,4, they find 

SI = b-'[Ss + 3c2] + 0{6^) , (5) 

5| = b-'^[S4 + 12c2S3 + ic3 + 124] + 0{5'^) , (6) 

where Sj are the galaxy amplitudes, Sj are the matter density moments, b = bi is the linear bias, 
and C2 = b2/b and C3 = b^/b are the lowest relative non-linear bias terms in equation (4). Since ^3 
and C2 may generally be of order unity, the contribution of the quadratic bias to ^3 is comparable to 
that from the linearly biased second-order skewness of the matter. Thus, it is inconsistent to assume 
either a purely linear bias or to ignore the gravitationally induced skewness, even in the limit of very 
small fluctuations. This method has been extended to scale-dependent bias models by Frieman and 
Gaztanaga (1994) and we make use of their results below. 

Detailed analysis of N-body results (Baugh, Gaztanaga & Efstathiou 1994 in preparation) show 
that the perturbation theory predictions for Sj agree remarkably well with the non-linear field for 
scales where ^2 — I5 while on smaller scales, ^3 in the simulations rises above the perturbative results. 
Thus, in Fig. 1, the comparison with perturbation theory must be restricted to scales larger than 
i? ~ 7 h~^ Mpc to be in the quasi-linear regime. This restriction also implies that we only need include 



Table 1: Best fit values of the non-linear bias parameters C2 and C3 for different values of b, fitted using 
the APM ^3 and ^4 and the predictions of different models. Each entry is a range of allowed values 
[Amax,Amin]x'^ whcrc x^ IS the gooducss of the fit at the ends of the range. 



r = o.2 


b= 1 


b = 1.5 


b = 2 


C2 
C3 


[-0.05,0.02]2.5 
[-0.3,2.0]6 


[0.42,0.49]2.5 
[-0.2,0.6]6 


[0.57,0.65]2.5 
[-0.0,0.8]6 


r = o.5 


b=l 


b = 1.5 


b = 2 


C2 
C3 


[0.13,0.20]4 
-0.2,1.8 6 


[0.29,0.38]2.5 
0.5,1.8]6 


[0.47,0.56]2.5 
0.6,1.9 6 


CGF 


b=l 


b = 1.5 


b = 2 


C2 
C3 


[0.44,0.48]ii 
[-0.3,1.0]i4 


[0.65,0.70]4 
[0.6,1.6]6 


[0.73,0.78]2.5 
[-0.0,2.1]6 



terms up to order J — 1 for Sj in the local bias expansion (4). Over this range of scales. Fig. 1 shows 
that the APM skewness and kurtosis are reasonably well fit (within the errors) by the low-density CDM 
model with no bias, i.e., bi = 1, 6jt>l = 0, whereas the unbiased CGF model does not agree with the 
data. 



We can use the five APM amplitudes S^iR) at R > 7 h~^ Mpc to find the best-fit quadratic bias 
parameter C2 for each model shown in Fig. 1, given a value of the linear bias b. Table 1 shows the range 
of values of C2 which give the best fit for each model. The goodness of the fit, parametrized by the 
X^(4) value, is shown as a subscript to the C2 interval, corresponding to 4 = 5 — 1 degrees of freedom. 
We then use a combined x test with the S^lR) and S^lR) APM values to constrain C3, given the range 
allowed for C2 in the S^lR) test for each model. In this case we have 9 data points at i? > 7 h^ Mpc 
(5 from ^3 and 4 from S4) and 7 = 9 — 2 degrees of freedom (since C3 is free and C2 is allowed to vary 
within a range). The best-fit ranges for C3 are also shown in Table 1 and correspond to a goodness of 
fit x'^i'^) — 6, except for the CGF model with 6 = 1 which has X^(7) c^ 14. 

For no linear bias, b = 1, the T = 0.2 model gives a good fit to the higher moments (at the 1-sigma 
level) with no non-linear bias, that is, for C2 = and C3 = 0, while the P = 0.5 standard CDM model 
requires C2 — 0.15 and C3 ~ 1, and even then only fits the data at the 2-sigma level. The CGF model 
does not fit the data with b = 1 for any choice of the c^. (x^(4) > 11). For & > 1, on the other hand, 
all three models can fit the data with different combinations of non-linear bias parameters C2 and C3 
as shown in Table 1. However, fits to the data with large values of b are in a sense ad hoc, because 
the agreement is obtained by suppressing the gravitational contribution to Sj and then fitting with the 
non-linear bias terms alone. These values do not really reflect agreement between the data and the 
model, but rather the possibility that, in any model, the observed, nearly constant moments can be 
reproduced by the constant non-linear bias terms in (5) and (6). 

In Fig. 2, we show the values of Sj obtained directly from the higher order correlations in the APM 
survey for J = 3 — 6. Again, there is reasonable agreement with the unbiased perturbation theory 
predictions for the P = 0.2 model and for the APM input spectrum from Baugh and Efstathiou (1993) 
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Figure 2: Hierarchical amplitudes Sj for J = 3 — 6 compared with the predictions in non-linear pertur- 
bation theory with no biasing. Curves correspond to the models in Fig. 1. 



on scales R > 7 h Mpc 



3 Conclusion 



The simplest interpretation of the higher order clustering data, in particular the agreement between the 
matter amplitudes Sj predicted in non-linear perturbation theory and the observed galaxy amplitudes 
Sj in the APM Survey (Figs. 1 and 2), is that there is no bias on scales R > 7 h~^ Mpc: the parameters 
in expansion (4) of 6g = f{6) are just 6 = 6i ~ 1 and bj^yi ~ 0. Otherwise, one either has to invoke a 
delicate, model-dependent balancing of terms in Eqs. (5) and (6) (as in the b = 1.5 models) or completely 
suppress the contribution of the matter amplitudes Sj via a large linear bias term b and generate the 
observed galaxy signal Sj with non-linear bias alone (e.g., for b = 2, the matter contribution to S^ is 
suppressed by a factor 4). As one goes to higher orders, this procedure must be done at each order. 

The predictions for the matter amplitudes Sj depend slightly on the shape of the matter power 
spectrum P{k). In Figs. 1 and 2 (solid curves) we have assumed that the galaxy spectrum measured in 
the APM survey is proportional to the matter spectrum. At large scales this spectral shape is similar to 
that of the T = 0.2 CDM model. For this case, Table 1 shows that 6 = 1 is compatible with no biasing, 
i.e. C2 = and C3 = 0, while b > 1 constrains the set of possible local non-linear bias parameters. 

One can also assume a different model for the matter spectrum P{k), such as the standard T = 0.5 
CDM model. In this case, it is necessary to introduce a non-local (scale-dependent) bias model, such 
as the cooperative galaxy formation model of Bower, etal., to explain the discrepancy in shape between 
the APM galaxy power spectrum and the matter density spectrum. We have shown here that scale- 



dependent bias alone, with no change in the overall amplitude (b = 1), is inconsistent with the observed 
higher order moments (Fig. 1 and Table 1). In principle, one could change the overall normalization of 
the spectrum (through a local bias) in this model; in this case, a reasonable agreement with the higher 
moments can be found for b c^ 2 hy fitting the values of C2 and C3. However, this choice of b disagrees 
with the COBE DMR normalization of the standard CDM model. 

Let us consider three models with initially Gaussian, scale-invariant fluctuations, to illustrate the 
imphcations of these results: standard CDM (SCDM) with Oq = 1, /i = 0.5; low-density CDM (LCDM) 
with 0.Q = 0.2, h = 1; and cold plus hot dark matter (CPHDM) with Qcold = 0.8, ^hot = 0.2, and 
h = 0.5. The linear matter spectrum P{k) in each model can be normalized using the COBE DMR 
measurements (Wright et al. 1994), which gives erg ~ 1 (SCDM, LCDM) and Gg ~ 0.7 (CPHDM). The 
APM normalization of the galaxy spectrum corresponds roughly to Cg ~ 1, so the linear bias for each 
model is fixed by COBE to bcoBE ^ 1 (SCDM, LCDM) and hcoBE ^ 1-5 (CPHDM). The shape of 
the CPHDM power spectrum is qualitatively similar to that of LCDM and to that inferred for APM 
galaxies. 

The COBE-normalized CGF model, based on SCDM and therefore with bcoBE — 1, is ruled out 
by the ^3 data in the APM survey, which gives X (4) > 11 (Table 1). As noted above, the LCDM 
model, also with bcoBE = 1, is compatible with no biasing, c^, = 0. The COBE-normalized CPHDM 
model, which corresponds approximately to the b = 1.5, P = 0.2 entry in Table 1, requires a particular 
combination of non-linear bias parameters to fit the ^3 and ^4 data, i.e., C2 — 0.45, C3 ~ 0.2. This 
judicious choice of non-linear bias parameters would need to be extended to C4 and C5, using the values 
of S5 and S(i in Figure 2. 

In the future, we expect that more accurate determinations of the large-scale higher order moments 
Sj and of the multi-point functions ^7(^1, ...,tj) (Fry 1994) will yield even tighter constraints on the 
nature of bias. 
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